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A new tool to determine masses and mass profiles using 
gravitational flexion 
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ABSTRACT 

In a recent publication, an aperture mass statistic for gravitational flexion was derived 
and shown to be effective, at least with simulated data, in detecting massive struc- 
tures and substructures within clusters of galaxies. Further, it was suggested that the 
radius at which the flexion aperture mass signal falls to zero might allow for esti- 
mation of the mass or density profile of the structures detected. In this paper, we 
more fully explore this possibility, considering the behaviour both of the peak signal 
and the zero-signal contours for two mass models-the singular isothermal sphere and 
Navarro-Frenk- White profiles-under varying aperture size, filter shape and mass con- 
centration parameter. We demonstrate the effectiveness of the flexion aperture mass 
statistic in discriminating between mass profiles and concentration parameters, and 
in providing an accurate estimate of the mass of the lens, to within a factor of 1.5 or 
better. In addition, we compare the aperture mass method to a direct nonparametric 
reconstruction of the convergence from flexion measurements. We demonstrate that 
the aperture mass technique is much better able to constrain the shape of the cen- 
tral density profile, obtains much finer angular resolution in reconstructions, and does 
not suffer from ambiguity in the normalisation of the signal, in contrast to the direct 
method. 
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1 MOTIVATION 

There is an ongoing debate regarding the exact shape of 
galaxy- and cluster-scale dark matter haloes in the Universe. 
N-body simulations carried out under the assumption of the 
standard concordance cosmological model, ACDM, suggest 
that cold dark matter haloes are well fit by Navarro-Frenk- 
White (NFW) density profiles (Navarro, Frenk & White, 
1997), which are defined by a characteristic virial radius, 
^200 (or, equivalently, the virial mass A/200), and a con- 
centration parameter, c. However, gravitational lensing ob- 
servations of field galaxies find that galaxy haloes consis- 
tently appear to show isothermal density profiles on the 
scales probed by strong lensing (Treu & Koopmans, 2002; 
Rusin, Kochanek & Keeton, 2003; Rusin and Kochanek, 
2005; Koopmans et al., 2006; Gavazzi et al., 2007; Czoske et 
al., 2008; Dye et al., 2008; Tu et al., 2009), resulting from 
the interplay between baryons and dark matter. There is fur- 
ther debate concerning the density profiles of galaxy clusters 
(Carlberg et al., 1997; van der Marel et al, 2000; Athreya 
et al, 2002; Katgert, Biviano & Mazure, 2004; Lin, Mohr & 
Stanford, 2004; Hansen et al., 2005; Lokas et al., 2006; Rines 
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& Diaferio, 2006; Wojtak et al., 2007; Okabe et al., 2009); 
although some gravitational lensing studies favour an NFW 
model over an isothermal model, there are claims that NFW 
models give a poor fit to some clusters (e.g. Newman et al., 
2009). 

In addition, there has been much discussion regarding 
the relationship between halo mass and NFW concentra- 
tion parameter. ACDM N-body simulations imply that the 
NFW profile is, in fact, a one-parameter (rather than two- 
parameter) model with the concentration parameter, c, be- 
ing related directly to the virial mass M200 as a power law 
with negative index (e.g. Navarro et al., 1997; Shaw et al., 
2006). Gravitational lensing observations involving galaxy- 
galaxy lensing in the Sloan Digital Sky Survey, for example, 
seem to support this supposed mass-concentration relation 
(Mandelbaum, Seljak, & Hirata, 2008). However, numerous 
gravitational lensing studies of clusters of galaxies, both in 
the strong and weak lensing regimes, have found significantly 
higher concentration parameters than would be expected 
from the mass-concentration relation (e.g. Abell 1689, for 
which a compilation of results can be found in Corless, King 
& Clowe, 2009). 

Furthermore, in a number of gravitational lensing stud- 
ies of clusters of galaxies, different sets of lensing measure- 
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ments are seen to disagree, and different (NFW) parametric 
models of the cluster with widely varying masses and con- 
centration parameters are seen by different authors to pro- 
vide good fits to the data. For example, gravitational lensing 
studies of the galaxy cluster Abell 1689 (see, e.g., Corless et 
al. 2009; Peng et al. 2009; and references therein) have found 
best-fit models that span a factor of 3.5 in mass and 3 in 
best-fit NFW concentration parameter. Such disagreement 
hinders efforts to compare cluster data with the predictions 
from N-body simulations described above, and thus to con- 
strain cosmological parameters. The technique presented in 
this paper is designed to assist in breaking some of these 
degeneracies, as even within a family of models such as the 
NFW model, the FM ap signatures for profiles with differ- 
ent lens masses and concentration parameters are seen to 
be substantially different from one another. 

Measurements of gravitational flexion, the second order 
gravitational lensing effects which give rise to skewness and 
arciness in galaxy images, have been shown by numerous 
authors to be adept at detecting galaxy- and galaxy group- 
size haloes both in the field and within clusters of galaxies 
(e.g. Okura, Umetsu & Futamase, 2007,2008; Leonard et 
al., 2008; Leonard, King & Wilkins, 2009), as well as pro- 
viding an alternative method by which the mass distribution 
within a cluster of galaxies might be constrained. In addi- 
tion, Lasky & Fluke (2009) have found that the first flexion 
signal from an SIS profile differs from that of an NFW sub- 
stantially at moderate separation between source and lens, 
whilst the second flexion signal shows strong variation when 
the NFW concentration parameter is varied. Flexion studies 
could therefore provide complementary constraints on halo 
profiles as well as their masses. 

Recently, Leonard et al. (2009; hereafter LKW09) devel- 
oped an aperture mass statistic for flexion, in direct analogy 
with that used for shear, and showed that it provided a ro- 
bust method by which structures within clusters of galaxies, 
on a range of physical scales, can be identified to appre- 
ciable signal to noise. This technique is formally identical 
to the standard shear aperture mass methods used in weak 
lensing (see, for example, Schneider 1996, Schneider et al. 
1998), but uses measurements of flexion rather than shear. 
This technique has the advantages that the noise properties 
of the aperture mass maps generated are very well under- 
stood, and that the filter functions used can be tuned to 
provide optimal signal-to-noise for a given lens profile. 

Moreover, the flexion aperture mass technique offers a 
robust method by which the mass distribution of a cluster 
of galaxies might be mapped out without the need for para- 
metric modelling. This means that the technique does not 
rely on assumptions about the mass profiles of the struc- 
tures responsible for the lensing signal, nor does it require a 
priori knowledge of the locations of the mass concentrations 
responsible for the lensing signal. It can therefore be used 
to place independent constraints on the mass distribution of 
clusters of galaxies, without the invocation of any assump- 
tions regarding the shape of the mass density profile of the 
structure. 

In addition, LKW09 noted that the zero-signal contours 
expected to be found around mass peaks vary in radius 
under change of mass profile, lens mass, filter shape and 
aperture size. LKW09 suggest the properties of this set of 
contours might thus be used to provide insights into the 



mass and profile shape of the structures identified using this 
method, thus offering the potential to discriminate between 
competing mass models in cases where degeneracies arise. 
In this paper, we more fully investigate this claim. 

We focus on two properties of the flexion aperture mass 
signal-namely, the peak signal associated with a given struc- 
ture and its zero-signal contour-and investigate their be- 
haviour for SIS and NFW profiles under variation in virial 
mass, NFW concentration parameter, aperture radius and 
filter shape in order to quantify the discriminating power 
of this technique. We demonstrate that the flexion aperture 
mass signatures of these two mass profiles are very differ- 
ent under changes to the aperture and filter properties. As 
a result of this divergent behaviour, aperture mass filter- 
ing of flexion data using a variety of aperture parameters 
provides a convenient and straightforward method for esti- 
mating both the mass and profile shape of the structures de- 
tected without the need for parametric modelling. Moreover, 
we demonstrate on simulated data that the total mass can 
be constrained to within a factor of ~ 1.5 for both galaxy- 
group and cluster scale haloes, and the input profile shape 
recovered, for a signal to noise in the flexion aperture mass 
measurement as low as 1. Finally, we demonstrate that this 
method significantly outperforms a Fourier transform-based 
direct inversion technique (similar to that used in Okura 
et al.'s 2008 analysis of Abell 1689), yielding both higher 
resolution and much better constraints on the shape of the 
cluster mass profile. 

This paper is structured as follows. In §[5] we provide a 
brief review of the origin of the flexion signal, an overview 
of the flexion aperture mass statistic, and a description of 
the mass models considered. In § [3] we consider the radial 
profile of the aperture mass for each of the mass models and 
for a range of filters and aperture radii, and for varying val- 
ues of the lens virial mass and concentration parameter (in 
the case of NFW lenses). The behaviour of the peak signal 
and the zero-signal contour under changes to the underlying 
mass model and changes to the aperture and filter proper- 
ties is considered. In § [4] we investigate the behaviour of a 
single mass profile under changes to the aperture and filter 
properties, and examine whether the use of several different 
combinations of aperture radius and filter shape in aperture 
mass reconstructions enables one to discriminate between 
different halo masses and density profiles. In § [5] we apply 
aperture mass and direct reconstruction techniques to flex- 
ion data obtained by raytracing simulations through a clus- 
ter extracted from the Millennium simulation (see LKW09 
and references therein) and compare the performance of the 
two methods. We conclude in § |S] with a discussion of our 
results and their implications for future flexion studies. 

Throughout the text, we assume a standard ACDM 
cosmology with Q m = 0.27, Qa = 0.73, and Ho = 
100/1 kms" 1 Mpc" 1 . 



2 BASIC FORMALISM 
2.1 Gravitational Flexion 

We begin by providing a very brief introduction to the flex- 
ion formalism. For a more complete discussion of the theory 
and description of the various techniques used to measure 
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flexion, the reader is referred to Goldberg & Bacon (2005), 
Bacon et al. (2006), Goldberg & Leonard (2007), Leonard 
et al. (2008) and Okura et al. (2007, 2008). 

Gravitational flexion arises when the lens field varies 
significantly over the scale of the lensed image. In this case, 
the lens equation - taken to be linear in weak lensing studies 
- must be extended to second order: 



fa ~ AijOj + -D 



ijkOjOk , 



(1) 



where f3 is the coordinate in the source plane and is the 
coordinate in the lens plane. These coordinates are related 
by (5 — — a(0), where a(0) is the deflection angle induced 
by the lens potential. A is the magnification matrix, which is 
related to the convergence, k, and the complex gravitational 
shear 7 = \j\e , and Dijk = dkAij is related to the first 
flexion, T — \T\e 1 ^ = <9k, and second flexion, Q — |(J|e 3 ^ = 
d-y, only. 

The first flexion signal gives rise to a skewness in the 
brightness distribution of the lensed galaxy image, and is 
a direct probe of the local gradient of the convergence. As 
such, it offers an ideal probe of cluster substructure, as it 
tends to be more sensitive to small-scale structures than to 
the large-scale cluster potential outside of the critical region 
of the cluster. 

The second flexion signal gives rise to a bending or arci- 
ness in the lensed galaxy image. Although we expect the sec- 
ond flexion signal to be larger for most mass profiles than the 
first flexion signal (see, e.g., Goldberg & Bacon, 2005; Lasky 
& Fluke, 2009), Goldberg & Leonard (2007) and Leonard et 
al. (2008) have found second flexion to be significantly more 
difficult to measure. We therefore restrict our discussion for 
the remainder of this paper to measurements of first flexion. 
The convergence is related to the first flexion through: 



J dV £>(x - x')J r (x) 



where 



and X — x\ + 1x2- 



E - = JF> 



+ ko , (2) 



(3) 



2.2 Direct Reconstruction Techniques 

To date, two groups have presented gravitational flexion 
studies of a cluster of galaxies with encouraging results. 
Leonard et al. (2008) presented the first measurements of 
flexion in a cluster of galaxies, and described a parametric 
modelling technique to reconstruct the distribution of mass 
within the galaxy cluster Abell 1689. Whilst this method 
showed good results in the periphery of the cluster, allowing 
mapping of the substructure content of the cluster, it was 
not able to constrain the central mass profile at all as a result 
of various masking techniques employed by the authors. 

Okura et al. (2008) demonstrated a nonparametric flex- 
ion reconstruction technique on the same cluster, with good 
results being obtained using only 5 background sources per 
square arcminute and a positive detection in the central re- 
gion of the cluster. Their method involves bin-averaging the 
flexion signal and reconstructing the convergence according 
to Equation(2] This is performed most simply and efficiently 



in Fourier space, where the relationship is expressed as (Ba- 
con et al., 2006): 



k\ + fe; 



(4) 



where X represents the Fourier transform of X. In taking 
the Fourier transform of k one recovers k — Ko, where ko is 
a constant of integration usually set by requiring the con- 
vergence to go to zero at large distances from the cluster 
centre. 

Such reconstructions may be problematic, however, as 
the resolution one can obtain in the output convergence map 
is strongly dependent on the background source density and 
distribution. Moreover, in the presence of significant sub- 
structure, bin-averaging over large pixels may have the effect 
of washing out the flexion signal. In addition, any masking 
schemes that result in a deficiency of sources in the central 
region of the cluster will result in the reconstructed conver- 
gence in this region being underestimated. It is for this rea- 
son that a parametric technique was chosen by Leonard et 
al. (2008) as a more appropriate method for reconstructing 
the mass distribution in Abell 1689 from flexion measure- 
ments. 



2.3 The Flexion Aperture Mass Statistic, F!M ap 

Aperture measures have long been used in weak lensing, and 
the effect of aperture filtering is to provide a robust measure 
of the mass distribution within lens systems whilst simulta- 
neously providing a straightforward method for determining 
the noise properties of the mass maps generated. In this way, 
aperture filtering offers an alternative nonparametric recon- 
struction method, which allows one to detect and map out 
mass concentrations within lens systems without requiring 
a priori knowledge or assumptions regarding the shape of 
the mass profile of the lens. 

The aperture mass statistic is defined as (Schneider 
1996; LKW09): 



i(x ) = J d 2 x k(x + xo) iw(|x|). 



(5) 



Using Equation [5] above, we can re-define this in terms of 
the measured flexion as 



m(xo) = 



j d 2 x w(x) J cfn' -Bjr(x — x' + xo)-F(x') 



+ Id x w{x)kq 



(6) 



This can be made independent of the constant quantity Ko if 
we require that the mass filter function w(x) is compensated, 
i.e. 



[ 

Jo 



x w(x)dx = 0. 



(7) 



Making the transformation y = x' — xo, and writing 
d 2 x = x dx d<j>, we obtain 

m(xo) = 
3? 



dyj(y + x ) / xw(x)dx / d<f> E* T {X - Y) 
Jo 



(8) 
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Applying the residue theorem to the integral over <j> yields 
(see Appendix B of LKW09 for a full derivation) 



m(xo) = — 5R 



d 2 y T(y - xo)tf 



'){x) dx 



"E-mode" (radially aligned) flexion as follows: 
m(xo) = j d 2 y JE(y;xo) Qriy), 



where 

Tb = 5ft [Je -10 ] 

= SR 



4 1 



(9) 



(10) 



(11) 



and 



--(») 



1 P 

- / a; tufas) dx. (12) 
y Jo 



We note that the E-mode flexion can equivalently be ex- 
pressed as 



Tn 



T-y 



(13) 



where here we have expressed the first flexion as a vector, 
rather than in complex notation, and y = y/|y|. 

As discussed in LKW09, any set of filters that are con- 
tinuous and satisfy Equations [7] and [T5] can be used in an 
FMa, p analysis. Ideally, one would choose a flexion filter func- 
tion, Qr(y), that traces the expected flexion signal in order 
to optimise the signal to noise in the FM ap maps generated. 
However, for the purposes of this paper, we assume that the 
mass profile is not known a priori. We choose the family of 
polynomial filters in LKW09, which are given by 



2(2 + or(l + 



1 - 



(14) 



tt! r(3 + o V R 2 , 

where 7? is the radius of the aperture being used. This family 
of filter functions has been shown by LKW09 to offer appre- 
ciable signal to noise for a range of different mass profiles, 
and is robust in the structures detected under variation of 
polynomial order I and aperture radius R. 

2.4 Mass Models 

We consider two of the most commonly-used models for 
gravitational lenses: the singular isothermal sphere (SIS) 
and Navarro- Frenk- White (NFW) profiles. Lasky and Fluke 
(2009) have presented an extensive description of the conver- 
gence, shear and flexion signals from both these models, as 
well as a Sersic profile, and we use their results - transformed 
into angular coordinates - in this paper. In addition, we as- 
sume our lenses to be circularly symmetric. Hawken and 
Bridle (2009) have presented a description of flexion from 
elliptical lenses; however it is sufficient for the purposes of 
this paper, and indeed makes the analysis somewhat simpler, 
to consider only lenses with circular symmetry. 

In order to compare these models, it is first helpful to 
define the virial radius, -R200, and the virial mass, M2oo- 
The virial radius is defined to be the radius within which 
the average density p(i?20o) = 200p c , where 

3H 2 



Pc = 



8ttG 



(15) 



is the critical density of the Universe, and 
3 f r 

= — p(x) x' 2 dx (16) 
r Jo 

for a spherically symmetric lens. The virial mass is then 
given by 



,, 8007T 3 

= — g — Pc-fl200- 



(17) 



We also remind the reader that the convergence, k, of 
a lens is related to its 2-dimensional projected surface mass 
density, E, by 



K{0) 



E(0) 



(18) 



where the angular coordinate, 0, is related to the physical 
projected position vector, £, and the observer-lens separa- 
tion, Dd, by — £/ 'Dd, and 



E c = 



c 2 D s 
4ttG D ds D d 



(19) 



where D s is the observer- source separation and Dds is the 
lens-source separatiorQ. 



2.4-1 Singular Isothermal Sphere Profile 

The SIS model has a 3-D density profile given by 



p{r) = 



2-wGr 2 



(20) 



where a is the velocity dispersion of the lens. The mean 
density within radius r is therefore given by 



p{r) = 



3a 2 



2nGr 2 

Setting this equal to 200p c , we find 

a 



which gives 



-R200 



M200 = 



50H 



2a" 



(21) 



(22) 



(23) 



/50GB 

Projecting onto the lens plane, the convergence is given by 



where Qe is the Einstein radius of the lens, given by 



(24) 



9 B = 4tt 



a \ Dds 



v v c J D B 

v c is the velocity of light, and the first flexion is given by 



(25) 



T{0) 



2I0I 2 



r0 



(26) 



All distances used are angular-diameter distances, unless oth- 
erwise specified. 
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2.4-2 Navarro-Frenk- White Profile 

The NFW profile has a 3-D density defined by 

ScPc 



p{r) 



(r/r s )(l + r/r s ) 2 



where 



200 



(27) 



(28) 



" 3 ln(l + c) - c/(l + c) ' 

and c is the concentration parameter, defined as c = 
R20o/r s . For a given choice of virial mass, the virial radius 
-R200 is computed from the definition in Equation 1171 

Before defining the convergence, it is useful to intro- 
duce some shorthand notation. First, we define x = £j/r s — 
\9\/0 s , where \6\ = £|/Dd- Second, we define a normalisa- 
tion factor: 

2p c S c r s 




(29) 



(30) 



(31) 



(32) 



In this notation, the convergence is given by 
where 



l(x) = 



and the first flexion is given by 

3 THE FM aD RADIAL PROFILE 



We now consider the expected behaviour of the FM &P signal 
for a given mass model, aperture size and filter shape. For 
convenience, we shall consider only apertures that lie along 
the a;-axis, i.e. the position vector from the centre of the lens 
to the centre of the aperture is given by xo = xoi. This is a 
valid simplification because the lenses under consideration 
all have circular symmetry, hence the radial profile of the 
FM ap signal along any given axis will provide a complete 
description of the expected signal. The displacement vector 
of a point within the aperture with respect to the centre of 
the aperture is denoted by y = y cos <f>i + y sin <f>j. Thus, the 
separation between the centre of the lens and a given point 
within the aperture may be expressed as 6 = xo + y = 
(xo + y cos (j>)i + y sin <j>j. 

The flexion vector always points along the direction of 
0, i.e.: 



H0) = \?\ iff 



(33) 



and, by definition, 

J"fi(y;x ) = ^— ^ 
V 



[xo cos 4> + y] . (34) 



s/xl +y 2 + 2x ycoscj> 
It is important to note that \T\ will, itself, generally be 



x ° 1 .0*10 4 




Figure 1. The expected -FM ap signal, evaluated numerically, for a 
lens of M200 = I0 15 h~ 1 Mq using a polynomial filter with I = 3 
and R = 120". The signal is plotted for an NFW profile with 
concentration parameter c = 3 and c = 12 and for an SIS profile. 
The dotted horizontal line shows m(xo)=0. 



dependent on xq, y, and (j>, as it is a function of 0. Thus, 
rewriting Equation 1 101 FM ap is given by 



i(x ) = 



2(2 + z) r (| + z) 



r(3 + /) 

\F\(x ,y,4>) [xp cos (j> + y] 
sj x\ + y 2 + 2x y cos <j) 



1-2- 
R 2 



1+1 



dy 



(35) 



As can be seen above, in general the flexion aperture 
mass is rather difficult to evaluate analytically, even for the 
simple case of a singular isothermal lens (unless the aperture 
is centred on the lens; this special case is discussed below). 
Therefore, in order to characterise its behaviour under dif- 
ferent mass models, aperture scales and filter shapes, it is 
necessary to evaluate the signal numerically. 

To do this, for each combination of virial mass, mass 
model, aperture size and filter shape, we consider 200 aper- 
tures evenly spaced in the range iro = [0",200"]. At each 
aperture location, we use a uniform random number gen- 
erator to generate 1000 y positions with —R ^ yi ^ R, 
and select those that lie within the aperture (i.e. those with 
y| $S R), so that the signal is evaluated on a uniform dis- 
tribution of points within the aperture. At each point, the 
E-mode flexion with respect to the centre of the aperture is 
computed, and the FMap signal at xo evaluated according 
to 



m(x ) = - V] T E .i Q(yi), 



(36) 



where n is the number density of points within the aperture. 

As only a finite number of points are being evaluated, 
this measure can be rather noisy. To reduce the noise, for 
each value of xo, the signal is estimated for 200 different 
sets of randomly-sampled data points, and the median value 
taken as an estimate of the signal. The resulting FM ap signal 
is plotted in Figure [1] for a sample of cluster mass lenses with 
M200 = 1O 15 /i _1 A/ , I = 3 and R = 120". 

It is immediately apparent that the FAf ap profiles differ 
greatly both in amplitude and in effective width. Therefore, 
comparison of these two features should provide information 
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about the underlying lens mass and density profile. We take 
the radius at which the signal goes to zero as a proxy for 
the effective width of the FM ap profile in each case. Other 
contours, such as the width at half the maximum signal or 
the radius of the minimum signal, could be used; however, 
we selected the zero-signal contour as both easily identifiable 
and convenient to fit. 

The profiles shown in the figure are not well-fit by a 
single analytic function. It is possible to fit these curves 
with a Chebyshev series of high order (~ 20; T. Nguyen, 
priv. comm.). However, the signal is equally well-fit (and, 
in some cases, better fit) by a piecewise continuous poly- 
nomial function, with a different polynomial behaviour seen 
for xo < xo.min and for xq > xo, mm, where xo >m in is the 
value of xo at which the minimum signal is seen. We find 
that the data in the region xo < xo, m in are generally well fit 
by a polynomial of 4 th order in xo- The zero-signal radius 
for an FM ap profile is therefore determined by finding the 
first real, positive root of the best-fit 4 th order polynomial. 

In order to assess the behaviour of these two fea- 
tures of the -FM ap signal for different masses, mass pro- 
files, aperture sizes and filter shapes, we carry out this nu- 
merical evaluation of the FM ap signal for W 11 ^ 1 M@ ^ 
M 200 < lO^ft^Me, R = [45", 60", 75", 90", 105", 120"], 
and I = [3, 4, 5, 6, 7, 8, 9, 10]. For each mass, we evaluate the 
signal for an SIS model, and for NFW models with concen- 
tration parameter c = [3, 4, 5, 6, 7, 8, 9, 10, 11, 12]. Note that 
throughout this paper, all angular scales are measured in 
arc seconds; therefore FM ap carries units of arcsec 2 . 

In addition, for all the simulations carried out here, we 
assume a lens redshift of Zd = 0.2 and a source redshift 
of 2 a = 1.0, typical of lens system geometries encountered 
in practice. These distances come into the normalisation of 
the flexion signals in different ways for the two mass pro- 
files considered here. The singular isothermal sphere model 
has a flexion measurement related to the source and lens 
separations through: 



•Fsis oc 0e oc 



-Dds 
D s 



(37) 



whilst the redshift dependence of the NFW flexion signal is 
given by: 



Jnfw oc 



PcTs 



H 2 {z d 



gdsgd 

D s 



(38) 



where H(zd) = HqWQmJI + Zd) 3 + Ha- These are simple 
scaling relations, however, and it is straightforward to com- 
pare measurements with different source and lens redshifts 
by appropriately scaling the FM ap signal. 



m(0) 



Noting that 



2(2 + Z)r(|+l) 

' r(3 + j) 
2(2 + r(|+Q 
; r(3 + Q 



l-V- 
R? 

r(2 + i) 
zr(|+0 



dy 



r(n + 1) = nT(n) , 
if n is a positive integer or half-integer, we find that 



m(0) = Be ( - + 1 ) R 



(40) 



(41) 



(42) 



Expressed in terms of the virial mass of the lens, the ex- 
pected peak FM ap signal is given by: 



m pe ak = m(0) = Aty 



5OGHM200Y D ds /5 , 
^1 ) "aU + 1 



(43) 

We now consider the point at which the FM ap sig- 
nal goes to zero. This manifests in flexion aperture mass 
reconstructions as a fairly well-defined zero-point contour 
(LKW09). It is not particularly straightforward to solve 
for this zero-signal radius analytically; therefore this radius 
is evaluated numerically by fitting the data in the region 
Xo < xo, m in with a 4 th order polynomial in Xo, and subse- 
quently finding the first real, positive root of this function. 

The zero-signal contour for the SIS model appears to 

have no dependence on mass. This is to be expected because 

2 

the aperture mass m(xo) oc 0e oc il/ 2 3 00 for all values of 
xo. Moreover, the SIS model has no mass-dependent scale 
radius. Therefore, in setting m(xo) = 0, the dependence on 
8e is removed, and one would expect the zero-point radius 
to be related only to the polynomial order I of the filter and 
the aperture radius R. The data imply that the zero-point 
radius is well-fit by 

Ro = — — y" (44) 

for all masses. 

Figure [5] shows a logarithmic plot of M200 vs 
fripeak/R (f + 1) for the entire numerical data set, consisting 
of 48 different combinations of I and R in each of 13 mass 
bins. The solid line shows the predicted signal from Equa- 
tion [43] The correlation seen here provides an excellent test 
of the numerical methods used to evaluate the FM ap signal. 
The figure also shows the behaviour of Ro vs R/(2l) 1 ^ 3 for 
the numerical data set. 



3.1 SIS model 

We now consider the behaviour of the peak FM ap signal 
and the zero-signal radius for the SIS mass profile under 
variation of the aperture size and filter polynomial order. 
The peak signal, by definition located at xo = 0, can be 
evaluated analytically for the SIS model. For this model, at 
xo = 0, we find 

Te = \T\(y) 



3.2 NFW Model 

The NFW profile is rather more complicated, making a fully 
analytic calculation of the peak FM a p signal for filters with 
polynomial order I > 0, or for a number of different values 
of I, R, and d s simultaneously, impractical and not partic- 
ularly instructive. Therefore, estimating the signal numeri- 
cally, and fitting the resulting data with a somewhat simpler 
function, seems a more appropriate approach. Indeed, this 
is necessary for the calculation of Ro, which cannot be com- 
puted analytically even for a simple lens profile such as the 
SIS. The peak signal and zero-signal radius are both found 
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Figure 2. Top Panel: The expected peak signal, normalised by 
aperture size R and filter polynomial order I, as a function of 
M200 for an SIS lens. The solid line shows the predicted signal 
from Equationl43l Bottom Panel: R vs R/(2l) 1 / 3 for an SIS lens. 
The solid line represents a 1:1 mapping. Note that for each mass 
bin plotted, 48 different combinations of I and R were used. Thus, 
the plot above represents 624 data points in total. 



Figure 3. Top Panel: The peak FM ap signal as a function of 
mass for an NFW lens with concentration parameter c = 3. The 
best fit curves of the form log 1 g(m pea fe) = a + 61og 10 (M2oo) + 
a!(log 10 (Af20o)) 2 are shown as solid lines. Bottom Panel: The zero- 
signal contour as a function of virial mass for an NFW lens with 
c = 3. The power law fits, Rq = AM^qq, are shown as dotted 
lines. 



to depend in a non-trivial way on the virial mass M200, 
concentration parameter (c), aperture radius R and filter 
polynomial order I. 

Figure [3] shows the behaviour of the peak signal and 
zero-signal radius as a function of M200 at fixed I and R for 
an NFW profile with a concentration parameter of c = 3. We 
find that the peak signal data are well fit by a polynomial of 
the form log 10 (m peaft ) = a + blog 10 (M 2 oo) + d(log 10 (M 2 oo)) 2 , 
whilst the zero-signal radius Ro appears to follow a power 
law of the form Ro = AM^oo- The figure also shows the best 
fit curves, the parameters of which are given in Table [T] 

Figure [3] shows the peak signal and zero-signal radius 
for fixed mass, concentration parameter (c = 3), and I as 
a function of R. Shown on a logarithmic scale, the be- 
haviour of the peak signal clearly appears to be a power 
law in R, i.e. m pea k = aR n . From the plots shown, it is 
apparent that a and n are strongly dependent on the mass 
and only weakly dependent on the polynomial order, I. Here 
again, the zero-signal radius appears to follow a power law 
behaviour, with Ro — bR p . Table [2] shows the fit param- 
eters for the combinations shown in the figure. The figure 
shows a degeneracy in values of Ro between the reconstruc- 
tions with I = 10, M200 = 10 15 /i _1 M and those with 



I 


R(") 


a 


b 


d 


l°gio( A ) 


n 


3 


45 


-9.6455 


1.3799 


-0.0356 


0.9846 


0.0355 


3 


120 


-12.0236 


1.7413 


-0.0456 


1.4476 


0.0293 


10 


45 


-7.7901 


1.1614 


-0.0286 


0.8093 


0.0383 


10 


120 


-10.4559 


1.5902 


-0.0416 


1.1886 


0.0381 



Table 1. Table of coefficients for the best fit curves for the data 
shown in Figure|3] the peak FM ap signal and zero-signal contours 
vs M200 f° r an NFW model with a concentration parameter of 
c = 3. The fits are of the form log 1 Q(m, pea fc) = a+61og 1( )(M20o) + 
d(log 10 (M 20 o)) 2 and R = AM£ Q0 . 



I = 3, M200 = 10 u ft l M . This is to be expected: for a 
fixed mass, concentration parameter and aperture size, in- 
creasing I will decrease Ro as the filter will be narrower in 
width. On the other hand, for a fixed concentration param- 
eter, increasing the mass increases the scale radius B s , which 
will have the effect of increasing Ro if all aperture parame- 
ters remain fixed. 

We now consider the behaviour of the peak signal and 
zero-signal contour as a function of I at fixed mass and R. 
Figure [5] shows this behaviour for various values of M200 and 
R, with c again fixed at 3. Again, this behaviour seems well 
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Figure 4. Peak FAf ap signal (top panel) and zero-signal radius 
(bottom panel) as a function of aperture size R for an NFW model 
with a concentration parameter of c = 3. The dotted lines show 
the best fit power laws. 



Figure 5. m pea k (top panel) and Ro (bottom panel) as a function 
of I for various choices of M200 and R for an NFW lens with a 
concentration parameter of c = 3. The dotted lines show the best 
fit power law to each set of data. The fit parameters are given in 
Table [3] 



/ 


M 2O o(M ) 


!ogio( a ) 


11 




P 


3 


10 11 


0.4311 


0.9576 


-0.1813 


0.9420 


10 


10 11 


0.4736 


1.1231 


-0.3007 


0.9227 


3 


10 15 


0.8408 


1.8893 


0.0041 


0.8927 


10 


10 15 


1.0908 


1.9300 


-0.1490 


0.9260 



Table 2. Table of coefficients for the power law fits to the data 
shown in Figure HI 'm pe(1 fe = aR n and Ro = bR p . Here, again, 
the signal is due to an NFW lens with a concentration parameter 
c = 3. 



fit by a power law in I with variable index; i.e. rn pea k = al n 
and Ro = bF . As expected, we see that increasing I while 
holding all other parameters fixed decreases the zero-signal 
radius. Table [3] shows the best fit parameters for this power 
law for the data shown in Figure 

Finally, we consider the behaviour of the peak FM ap 
signal and the zero-signal radius for an NFW profile of fixed 
mass, with fixed aperture parameters R and I, as a function 
of concentration parameter c. The data have a strong depen- 
dence on halo mass, therefore Figure [6] shows the behaviour 
of the peak signal and zero-signal contours for various com- 
binations of I and R for both a halo of mass lO 11 /!" 1 M Q and 
10 15 h _1 M Q . Here, again, the behaviour of each appears to 
be a power law of the form m pe ak = ac n and Ro = be? . 
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Table 3. Best fit parameters for the data shown in Figure [5] 
assuming the relationships to be of the form m pea j. = al n and 
Ro = bl p . The lens in this case is an NFW lens with a concentra- 
tion parameter of c = 3. 



Table 3] shows the best fit parameters for these power laws 
for the data shown in Figure [6] 

The behaviour of the peak FM ap signal from the NFW 
profile differs greatly from that of the SIS model as the aper- 
ture radius and filter polynomial order are changed. More- 
over, the behaviour changes with halo mass and concen- 
tration parameter in a non-trivial way. A similar trend is 
seen when considering the behaviour of the zero-signal ra- 
dius, where the slope of the power law behaviour changes 
based on aperture size, filter shape, virial mass and concen- 
tration parameter. Again, it is not possible to model this 
behaviour simply and arrive at a general analytic expres- 
sion for the zero-signal radius that simultaneously describes 
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Figure 6. The figure shows the behaviour of the peak FM ap signal (top row) and the zero-signal radius (bottom row) as a function of 
the NFW concentration parameter for fixed values of I and R for a low mass halo (left panel, 10 11 ?! - 1 Mq), and a high mass halo (right 
panel, lO 15 ^" 1 Mq). The dotted lines show the best fit power law to the data, the fit parameters of which can be found in Table l4l 



M 20 o(M ) 


I 


R(") 


lo Sio( a ) 


n 


l°Sio( fe ) 


P 


10 11 


3 


45 


1.1171 


0.1986 


1.3931 


-0.0435 


10 11 


3 


120 


1.6694 


-0.0980 


1.7821 


-0.0149 


10 11 


10 


45 


1.3612 


0.3521 


1.2545 


-0.0693 


10 11 


10 


120 


1.9908 


0.0222 


1.6273 


-0.0328 


10 15 


3 


45 


2.4788 


1.1720 


1.5705 


-0.1092 


10 15 


3 


120 


3.3976 


0.9483 


1.9467 


-0.1098 


10 15 


10 


45 


2.6293 


1.2265 


1.4302 


-0.1013 


10 15 


10 


120 


3.5409 


1.0588 


1.8376 


-0.1427 



Table 4. Best fit parameters for the data shown in Figure [6] assuming relationships of the form m pe 



and Rq = bc p . 



its behaviour as a function of all the parameters one might 
vary. 

This vastly differing behaviour implies that if one were 
to consider several flexion aperture mass reconstructions of a 
lens field using different aperture radii and filter polynomial 
order, one might be able to distinguish not only between an 
SIS model and an NFW model, but between NFW models 
with different masses and concentration parameters. More- 
over, considering the behaviour exhibited in Figures [4] and [5l 
it appears that a change in aperture radius shows a greater 
change in the overall peak signal than a change in polyno- 
mial order, and thus might provide a better discriminant 
between profiles than the filter polynomial order. 



4 MODEL-SELECTION WITH M ap 

Having characterised the expected behaviour of the peak sig- 
nal and zero-signal radius of the FM ap -statistic for the SIS 
model and a family of NFW models of varying concentra- 
tion parameters, one might well ask how the differences seen 
between the various models might be used to constrain lens 
parameters in real observations. The first key feature seen in 
§|3]is that the value of the peak signal is a strong function of 
the virial mass of the lens, regardless of which mass model 
one is considering. This implies that a measurement of the 
peak FM ap signal will automatically, at very least to within 
an order of magnitude, give an indication of the mass of the 
structure responsible for the lensing signal. 
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Further, the data above demonstrate that, under a 
change in aperture radius, the peak signal and zero-signal 
radius change in different ways for each of the two models, 
and a strong variation is also seen between NFW models 
with different concentration parameters. Such divergent be- 
haviour is also seen when varying the filter polynomial order, 
I, though to a lesser extent. The simulated data therefore 
suggest that given a set of FM ap data for various combina- 
tions of I and R, the mass and mass profile of a halo might 
be determined uniquely by considering the behaviour of the 
peak signal and zero-signal radius as I and R are varied. In 
other words, -FM ap filtering of the flexion signal for a range 
of aperture parameters allows one to delineate between var- 
ious mass models, offering a method by which degeneracies 
between parametric models might be broken. 

In order to assess this further, we utilise the simu- 
lated radial profiles described in § [3] above. These simula- 
tions were evaluated for M2oo/(10 11 /! -1 M ) = 1, 2, and 
5 xl0°, 10 1 , 10 2 and 10 3 . For each mass, an SIS pro- 
file was evaluated, as well as NFW profiles with c sam- 
pled in integer values in the range 3 c ^ 12. Further, 
the -FM ap signal for each model was evaluated using all 
possible combinations of / = [3, 4, 5, 6, 7, 8, 9, 10] and R = 
[45", 60", 75", 90", 105", 120"]. In order to improve the reso- 
lution attainable at the low-mass end, additional simulations 
were carried out for the SIS and all the above NFW profiles 
with A/200 = [3,4,6,7,8,9] x lO 11 /^ 1 M Q and all combi- 
nations of I = [3,5,7,10] and R = [60", 90", 120"]. Simi- 
larly, at the high-mass end, we added simulations involving 
M200 = [3,4,6,7,8,9,20] x 1O 14 /i _1 M and all combina- 
tions of I = [3, 5, 7, 10] and R = [45", 60", 75", 90", 120"]. 

The discriminating power of the FA/ ap -statistic is tested 
as follows. We consider a lens of a given mass M200 and con- 
centration parameter c*. We assume that we are able to 
measure the FM ap -statistic for a number of combinations 
[I, R] with a signal to noise of S and an error on the mea- 
surement of Ro given by or . For each combination of / and 
R, we find the expected value of m* eak and Ro for our test 
model from the simulations already carried out. We then 
scan the parameter space to find any other models with val- 
ues of m pea h and Ro that fall within our error bars for the 
combination of I and R being used; i.e. m peak — (J pea k ^ 
mp ea k < rn* eak + <j pea k and Rf, — or < Ro < Ro + or 
(note that a vea k = Tn Pea k/S)- As I and R are changed, the 
lists of models that fit the test model's "data" within the 
error bars are compared, and only those models found to be 
compatible with all the "measurements" carried out using 
our model lens are retained. 

4.1 Galaxy-Mass Halo 

As a first test of this method, we consider a theoretical lens 
with a mass of M 2 *oo = 5 x lO 11 /!" 1 M with an NFW pro- 
file with concentration parameter c* = 12. We assume that 
we are able to measure the FM ap -statistic for I — [3, 5, 7, 10] 
and 7? = [60", 90", 120"] with a modest signal to noise of 2.0 
and an error on the measurement of Ro of ±1". After iterat- 
ing over all the available combinations of I and R, the only 
models in the simulated data set that are found to have con- 
sistently fallen within the error bars on the theoretical mea- 
surements are: M 20 o = 4x lO 11 /). -1 M with c = 7, 8, or 12, 
M200 = 5 x 10 u ft -1 M s with c = 10, 11, or 12, and 



M200 = 7 x lO 11 /! -1 A/ Q with c = 12. If the error on the 
measurement of Ro is reduced to ±0.5", the degeneracy is 
completely broken and we recover our input model, even for 
a signal to noise of 1. Note that these results are typical for 
this mass range. Changing the input value of c*, or using an 
SIS profile, showed a similar success rate. 

We can conclude therefore that this method is successful 
at recovering both the total mass and input concentration 
parameter to within a factor of < 1.5, and effective at elimi- 
nating alternative parametrisations for the mass profile such 
as the SIS model or NFW models with significantly higher 
or lower concentration parameters, given a sufficient num- 
ber of FM ap realisations. Moreover, in the low-mass case, 
reduction in the errors on measurements of Ro are most ef- 
fective at reducing the number of models compatible with 
the data, and thus the error on the estimates of the mass 
and concentration parameter. 

4.2 Cluster-Mass Halo 

We now consider a lens with a mass M 2 * 00 = 5 x 10 14 /i _1 M Q 
and concentration parameter c* = 4, and consider mea- 
surements at I = 3, R = [60", 90", 120"] with S = 2 and 
gro = 1". After only three iterations, the field of possible 
models has been narrowed down to M200 ~ 2 x 10 14 /i -1 M 
with c = 3, M200 = [4, 5, 6] x 10 14 ft -1 M Q with c = 4 and 
A/200 = 7 x 10 14 /i _1 Mq with c = 5. Increasing the signal 
to noise to 2.5 isolates models with c = 4 and A/200 in the 
range 4 x 10 14 /i _1 A/ ^ A/ 20 o < 6 x 10 14 /i _1 A/ . In this 
case, reducing the error on Ro to 0.5" does not completely 
break the degeneracy, but does reduce the number of avail- 
able models to 2, finding A/ 20 o = 4 - 5 x 10 14 /i _1 M Q with 
C = 4. However, increasing the number of models iterated 
over does break the degeneracy for modest error estimates. 
Using I = 3 and R = [45", 60", 90", 120"] with 5 = 2.5 and 
oro — 1.0 returns our original model. 

In the case of high mass haloes, far fewer FA/ ap real- 
isations are required for convergence than in the low mass 
case. The data here suggest that with only four FM ap re- 
constructions of a cluster-mass halo and reasonable errors 
on the zero-signal radius and the peak signal to noise allow 
us to recover our input model uniquely up to the resolu- 
tion of our sample data set. In this case, however, greater 
benefit is obtained by increasing the peak signal to noise in 
the FMap reconstruction or increasing the number of recon- 
structions obtained, rather than decreasing the errors on the 
measurement of Rq. 



5 TEST ON A SIMULATED CLUSTER 

In order to more realistically evaluate the performance 
of this method, we employ the i*M ap data from LKW09 
for a simulated cluster taken from the Millennium simu- 
lation (Springel et al., 2004). This cluster, on extraction 
from the simulation, was found to have a virial mass of 
M 2 * 00 = 1.23 x lO 15 ft _1 A/ , and is located at a redshift 
of 2d = 0.21. Background sources were generated at ran- 
dom positions with a density of 35 arcmin~ 2 at a redshift of 
z s — 1.0. The background sources were assumed to have zero 
intrinsic flexion, and an ellipticity drawn randomly from a 
Gaussian distribution with a standard deviation of a t = 0.2. 
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Figure 7. The fraction of pixels containing at least 7V m i n sources 
as a function of pixel size. 



These galaxies were lensed through the input cluster 
potential using a ray-tracing procedure (see LKW09 for de- 
tails of this calculation). As the galaxies had zero intrinsic 
flexion, the dominant sources of noise in the reconstructions 
arise from measurement errors and errors in the ray-tracing 
code resulting from the high degree of pixellation of the in- 
put convergence, which was not smoothed. As a result, one 
would expect a somewhat lower signal to noise to be seen in 
real data assuming the same density of background sources. 



5.1 Direct Nonparametric Mass Reconstruction 

When carrying out a direct reconstruction, the first issue to 
be addressed is that of the appropriate choice of bin size over 
which to average the flexion signal. A balance must be struck 
between using a fine grid, in order to resolve structures on 
both small and large scales, and having a statistically sig- 
nificant number of background sources within each bin. In 
Figure we show the fraction of pixels containing at least 
A^min sources as a function of bin width for N m i n = 1 and 3. 

To alleviate additional pixel noise in our reconstruc- 
tion resulting from incomplete coverage, we choose a bin 
size of 0.43 arcmin per pixel, resulting in a reconstructed 
convergence map of 16 x 16 pixels. In order to reduce edge 
effects when carrying out the Fourier transforms, we add 
zero-padding on all four sides of the image such that the 
padded image consisted of 128 x 128 pixels. To estimate the 
noise in our reconstruction, the flexion vectors were rotated 
by a random angle, and the reconstruction repeated. This 
procedure was carried out 1000 times, and the error at each 
pixel location estimated by the standard deviation Oi of the 
randomised reconstructions. 

Finally, the reconstructed convergence must be ap- 
propriately rescaled. As our data covers only the central 
lh^ 1 Mpc of the cluster, we do not expect the convergence 
to go to zero at the edge of the cluster. We therefore set the 
integration constant ko such that the weighted mean of the 
reconstructed convergence, defined as: 
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Figure 8. Top Panel: A reconstruction of the convergence map 
of the cluster from flexion measurements. Bottom Panel: A com- 
parison of the radially-averaged convergence profile of the input 
convergence (dotted line with points) and the reconstructed con- 
vergence (open diamonds with error bars). 



is equal to the mean convergence in the input map, found 
to be k = 0.13. 

Figure [8] shows the convergence reconstructed in this 
manner, as well as a comparison of the radial profile of the 
reconstruction with that of the input convergence map. Note 
that the error bars represent the 1-a errors on k — ko • Mod- 
elling the convergence as a power law of the form k oc r~ a , 
we find that the convergence reconstruction shows a shallow 
slope (a = 0.19), while the input map shows a much steeper 
slope (q ~ 0.66). We note also that the input convergence 
appears to follow a broken power law, rather than having 
a constant power law index, whilst the reconstruction does 
not appear to follow this behaviour. 

The underestimation of the central slope is not entirely 
unexpected. The reason is that one expects a decrease in the 
number of background sources with flexion measurements 
in the central regions of the cluster, resulting from sources 
being lensed away from the centre, and from blending of 
background images with those of cluster members. Figure [9] 
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Figure 9. The background sources averaged in annuli of width 
0.5 arcminutes as a function of distance from the centre of the 
cluster. 



shows the mean background source count evaluated in annuli 
as a function of distance from the cluster centre, demonstrat- 
ing a significant underdensity in the central 1.5 arcminutes 
as compared with the periphery. The lack of sources in the 
central region limits the ability of this method to constrain 
the potential in that regime. 

Note that due to our choice of normalisation constant, 
the integrated mass found through the direct reconstruction 
method will be consistent with that of the input map; how- 
ever, it is clear from Figure[5]that this method will not allow 
one to determine with any confidence the true shape of the 
underlying mass profile. Moreover, the choice of normalisa- 
tion constant relies either on the presence of data at large 
radii from the centre of the cluster, or on other external in- 
formation. In the absence of either of these, the choice of 
normalisation, and hence the cluster mass thus derived, is 
somewhat arbitrary. 



5.2 FM ap Analysis and Results 

As in LKW09, the _FM ap -statistic was measured for I = 5 
and R = [60", 90", 120"]. Hence, to test how well this profile 
might be identified, FM & p data for an SIS model and the full 
range of NFW models were obtained for MJoo, an d added 
to our sample of test profiles. 

The peak signal to noise in the reconstructions was 
found to be S = 2.5, 3.0, and 3.3, in order of increasing 
aperture radius, and we assume an error on the zero-signal 
radius measurements of <tr — 10", which corresponds to 
~ 5 pixels on the input convergence map. 

We note that in each reconstruction shown in LKW09, 
there is a significant peak seen at the top edge of the field of 
view. In the R = 60" FM^p reconstruction, this feature actu- 
ally obtains higher signal to noise than the central peak. This 
feature is believed to arise due to edge effects, as the signifi- 
cance of this peak is reduced in subsequent reconstructions, 
and its position changes somewhat as the filter is changed. 
This indicates a noise feature, rather than a true signal. As 
we are considering a radially-averaged signal here centred on 
the centre of the FM ap field, and this particular feature is 
located far outside the zero-signal radius of the central mass 



peak, this feature will not contaminate our measurements 
substantially. 

Using the three FM ap reconstructions of the N-body 
cluster, the data indicate that the radial profile of the clus- 
ter is well-fit by an NFW mass profile with c = 3 and 
M 200 = 1, 1-23 or 2 x 10 15 /i _1 M Q . The radial profiles for 
the FM ap signal from the cluster for each aperture radius 
are plotted in Figure 1101 along with the radial profiles of 
the three best-fitting models. The effect of the cluster's el- 
lipticity and associated substructure can be clearly seen in 
its radial profile, with a secondary peak being seen around 
xo = 130" in each reconstruction. 

All three models show values of R® that are smaller than 
that seen in the cluster data, implying that the best fit model 
might be one with an even smaller concentration parameter 
(see Figure HJJ. However, the peak signal does appear to be 
best fit across the three reconstructions by the profile with 
M200 = 1-23 x 10 15 /i _1 M Q , as expected. 

Furthermore, we can compare this best-fit models to 
the radial profile of the convergence map. Figure [TT] shows 
the measured radial profile of the convergence of the lens. 
The dotted curves show the convergence profile for the three 
best-fitting mass models described above. Clearly the two 
lower-mass models offer excellent fits to the cluster data, 
thus demonstrating the power of this method to accurately 
characterise the mass profile of a lens using a relatively small 



6 DISCUSSION 

In this paper, we have presented the expected FM^p sig- 
nal from two common mass profiles: the Singular Isothermal 
Sphere model and the Navarro-Frenk- White model. This sig- 
nal is characterised by a peak value, found when the aper- 
ture location coincides with a mass concentration, and a 
zero-signal contour. By simulating radial -FM ap profiles for 
a large number of combinations of halo mass, NFW concen- 
tration parameter, aperture radius and filter polynomial or- 
der, the behaviour of the peak signal and zero-signal contour 
has been characterised for both the SIS and NFW density 
profiles under the assumption of circular symmetry. 

The behaviour of both these measures is seen to diverge 
rather sharply between different models with the same mass 
as the aperture properties are varied. This implies that one 
might be able to use the FMap-statistic to constrain the 
masses and mass profiles of structures detected simply by 
varying the aperture parameters and noting the change in 
the peak signal and the location of the zero-signal contour. 

Indeed, it has been shown using simulated data from an- 
alytic models that with modest signal to noise and relatively 
large errors on the measurement of the zero-signal radius, it 
is possible to isolate an input model using a reasonable num- 
ber of combinations of aperture size and filter polynomial or- 
der. It also becomes easier to discriminate between masses 
at the high-mass end of the spectrum considered here, with 
higher mass haloes requiring fewer combinations of I and R 
to reach convergence. 

This was further reinforced by testing the method on 
FMap data using the simulated cluster of LKW09 as the in- 
put lens. This lens has a virial mass of 1.23 x 10 15 /i _1 Mq, 
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Figure 10. The figure shows the radially-averaged -FM ap signal 
for the N-body cluster described in LKW09 (solid curves), and 
the simulated data corresponding to the two best fit models from 
the simulations carried out in this paper (discrete points) for / = 5 
and R = 60" (top panel), 90" (middle panel) and 120" (bottom 
panel). 

was able to recover this mass to within a factor of < 1.5, and 
to correctly identify the shape of the profile as an NFW with 
low concentration parameter. This success, despite the cir- 
cularisation of a clearly elliptical lens profile, demonstrates 
solidly the power of the _FM ap -statistic to discriminate be- 
tween haloes of different masses. 

Moreover, this method shows a dramatic improvement 
over direct nonparametric reconstruction techniques. Firstly, 
the resolution that one can obtain using a direct method 
is limited to a far greater extent than aperture methods. 



Figure 11. The circularly averaged profile of the convergence of 
the N-body lens described in the text. The dotted curves shows 
the expected profile for an NFW lens with c = 3 for various lens 
masses as illustrated in the legend. 

This is because the number of bins one can use is limited 
by the density and distribution of background sources. We 
find, for the simulated cluster, a best resolution of just under 
half an arcminute. This coarse resolution is not problematic 
with the simulated cluster presented here, as it has a fairly 
smooth large-scale distribution. However, where there are 
significant substructures (such as in Abell 1689), averaging 
a signal over large bins has the effect of washing out the 
signal from smaller structures. 

Aperture measures, on the other hand, filter the signal 
through a large aperture (R > 0.5 arcminutes, in general), 
but are able to obtain fine resolution in the output aperture 
mass by spacing the apertures close together. Thus, finer 
structure can be resolved using these methods. 

In addition, we have shown that the ability of direct re- 
construction techniques to accurately characterise the shape 
of the central density profile in clusters is very limited, re- 
sulting from a decrease in the number of background sources 
with flexion measurements in the central regions of the clus- 
ter. Since direct methods rely on the presence of reliable 
flexion measurements within each pixel, the ability to mea- 
sure the central density profile is generally rather limited. 
Again, as aperture measures involve filtering the signal from 
a large area, they are less limited by a dearth of sources in 
the central regions. 

Lastly, whilst the direct method needs to be adjusted 
by a somewhat arbitrary constant shift, the aperture mass 
statistic is, by construction, independent of this normalisa- 
tion. 

There are often found to be somewhat large differences 
seen in mass estimates of clusters of galaxies found using 
strong lensing, weak lensing, or a combination of the two. 
The flexion aperture mass statistic offers a method to break 
the degeneracy between these models, as it clearly discrim- 
inates between different masses and mass models. As a re- 
sult of filtering the data, different mass models exhibit sig- 
nificantly different behaviours. Thus, this method offers a 
unique and novel model-selection method that is comple- 
mentary to other techniques currently used, and requires no 
prior assumptions about the underlying mass distribution. 

The issue of finding the best-fit mass profile is important 
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to cosmology. Cosmological simulations suggest dark matter 
haloes should follow a universal (NFW) density profile with 
concentration parameter anti-correlated with mass. There 
is much discussion regarding whether observations support 
these predictions. Thus, the ability to break degeneracies 
between models is important to cosmology; in this paper, 
we offer a clear prescription for a method by which this can 
be achieved. 

There is much work remaining to be done, of course. 
A more finely-sampled parameter space would allow for in- 
creased accuracy in model selection, improving the signal 
to noise in the measurements (by improving the filter func- 
tions used, for example) would increase the discriminating 
power of the _FAf ap -statistic, and considering the expected 
radial signal from elliptical lenses or models involving sub- 
stantial substructure would provide a more realistic sample 
of model templates from which to select. However, the work 
presented in this paper represents an excellent first step, 
and highlights the promise that the aperture mass statis- 
tic from flexion shows to be able to discriminate between 
masses, mass profiles and NFW concentration parameters 
using only a few combinations of aperture parameters. 
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